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Abstract 

A formalism is presented to obtain closed evolution equations for asymp- 
totic probability distribution functions of turbulence magnitudes. The for- 
malism is derived for a generic evolution equation, so that the final result 
can be easily applied to rather general problems. Although the approxima- 
tion involved cannot be ascertained a priori, we show that application of the 
formalism to well known problems gives the correct results. 
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I. Introduction 



Probability distribution functions (pdf's) are powerful tools for the ex- 
perimental and theoretical study of turbulence. Experimental pdf's can be 
obtained rather directly from the measurements of different magnitudes, and 
are at the same time a very convenient tool for the representation of such 
results. From the theoretical point of view, the equations for the evolution of 
pdf's can be readily obtained from the basic dynamical equations, although 
generally not in closed form. Closure requires the determination of condi- 
tional averages in terms of the pdf itself, which sometimes can be done based 
on additional information given by experiments,^ and numerical simula- 
tions.^ Closures of this kind have allowed to derive very important results 
directly from the original dynamical equations, such as the limiting and time 
dependent pdf of passive^ and non-passive^ scalars. We present here a pro- 
cedure for closure of the pdf equation, at least of its asymptotic form, that 
can be applied in a systematic way to general systems. It is important to 
mention that one can only show that the model is applicable if a reasonable 
conjecture is valid. We show that the model gives correct results in vari- 
ous systems whose behavior is known. The advantage is that it is very easy 
to apply and so it is valuable to attack new problems for which no useful 
information exists. 

II. Formalism 

The basic idea is most easily introduced considering a generic scalar field 
4>(x,t), where t is the time, and x denotes the spatial coordinates. The field 
satisfies a generic evolution equation of the form 



in which iV [(f)] is a spatial functional of (f> t) evaluated at the actual time 
t, and / is a stochastic forcing introduced to allow eventual consideration 
of statistically stationary pdf's. N [(f)] contains in general linear and non- 
linear terms, and we restrict the non-linear terms to entire powers of (f>, not 
necessarily local. The induced probability distribution of a given functional 
of 0, G [(f>] , can be conveniently represented as 




(1) 



?(0 = (i(GM-0), 



(2) 
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in which (...) means average over ensembles of realizations, and 5 is Dirac's 
delta function. The evolution equation of P (£) is then determined from Eq. 
© as 



OP d 
~dt ~ ~d£ 



,H; *[N[<j>,x 1 ]+f(x 1 ,t)} 



(3) 



5(f) (xi) 

where ( ...| £) represents average conditioned on G [0] = £, and the asteric in- 
dicates integration over repeated variables. We will consider only functionals 
G [0] that are linear in 0, so that (5G/5<p) * will be an operator independent 
of which we denote as 

7T7 x = h {x x ) . 
b(p{xi) 

Besides, as h (xi) is independent of 0, only (N [0] + f\ needs to be evalu- 
ated. 

To evaluate (N [<f>] \ £) let us consider the simplest linear term, 

T 1 (y,0 = {<l>(y)\h(xi)*<l>(x 1 )=0 

where y is an arbitrary spatial position and, from now on, no explicit indi- 
cation of the time is made. The exact, formal expression of Ti(y,£) can be 
written as 

r( n = /^WW(MsiW(si)-flP^ U) 

jT[<f>]5(h(x 1 )*<f>(x l )-0T><f> ■ 1 ' 

Here Dirac's delta function S is used to select only those fields that satisfy the 
condition h(x{)* <f>(xx) = ^P[4>\ is the probability density functional (pdF) 
of the field 0, and the integrations are meant to be functional integrations 
over fields 0, with an appropriate measure Z>0. In principle, given the system 
(JTj) a Martin-Siggia-Rose lagrangian can be determined in terms of which to 
express J-[(fi\. Once this is done, solutions of (j4]) can be obtained in terms 
of perturbative expansions, or non-perturbative approaches such as rapid 
descent type of integrations around significant field configurations. We follow 
here the perturbative approach, setting up an infinite series that represents 
(H]) in a formally exact manner, and then dividing this series into two infinite 
other series, such that one of them is term by term much larger than the 
other in the limit of large £. We keep only these larger terms and sum the 
resulting infinite series. 

To proceed we write both integrands appearing in (j2J) as series involving 
Gaussian functionals. For this we write 

^[0]= Jo[0]£[0], (5) 
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where 

1 



^,[0] = 1 exp[-/t(xi) * 0(xi) - -cr(xi, ar 2 ) * 0(xi)0(x 2 )], (6) 

is a Gaussian functional, with k(x) and cr(x, a;') functions to be determined, 
and Zq a normalization constant ensuring that J J r o[4>] T^'P — !• is 
in general a non-Gaussian functional factor defined by the very expression 
(j3J), satisfying also the normalization condition J •7 r [</>] T><p = 1. We further 
expand £?[</>] clS cl functional Taylor series 

£[0] = "T^ Xl ' X «) * 0(^l)---0(^n), (7) 



n! 

n=0 



where the n = term is to be understood as constant. 

We turn now to the evaluation of (j3j). For this, note that the denominator 
of this expression is the pdf (J2J), P(£), which, using the Fourier representation 
of Dirac's delta function, can be written as 



p(0 = yVM<KM*iW(zi)-£W 

(27T)" 1 / dkexp{-ik£) 



J r [<t>]G[<t>) exp [ifc/i (a*) * 0(a?i)] D0. (8) 
It is now useful to introduce an auxiliary functional given by 
A[J] = (27r) _1 J dkexp(-ik£) 

JF o [0] exp [J (xi) * 4>{xi) + ikh (xi) * <f>(xi)] T>(f>. (9) 



A plays the role of a potential from which (jSJ) and the numerator of (jl]) can 
be obtained by differentiation as it is immediately seen from (Q that, for 
instance, 

SA 

«H = (10) 

and so, successive derivatives allow to produce the factors appearing in the 
series of Q[<f>] (see (j^)). In this way, we write the denominator and numerator 
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of fll]) respectively as 

p(0 = J £^ K ^-> x ^* A{n) \j=o> ( n ) 



n! 



and 



where 



ri( ?/> e)p(e)=y;A^i,-^n)* 



it 



(12) 



j=o 



" 8J{x l )...8J{x n )' 1 ; 

Evaluating the Gaussian integrals in (jSJ) and ([9]) we obtain 

A[J] = iVexp (-A~ x £ 2 /2) exp W[J - k], (14) 

where iV is a factor independent of J and £, A = £(xi, 22) * h(xi)h(x2), and 

= 1/2 x 2 ) * ^(xi)X(x 2 ) - A^b 2 + 2A~ 1 6£] , (15) 

with b[X] = S(xi,X2) * h(xi)X(x2) a linear functional of X, and S is the 
inverse of a in the sense that 

S(xi, x 2 ) * cr(x 2 , x 3 ) = cr(xi, x 2 ) * S(x 2 , x 3 ) = <5(xi - x 3 ). (16) 

In this way, (JI]) can be calculated by explicit differentiation through ( ITTjl and 
( fl2|) . the series so obtained being formally equivalent to the original integrals. 

So far we have only set up a generic series expansion in analogy with 
usual methods, and to obtain useful results this series needs to be evaluated 
at least approximately. We assume no intrinsic small parameter, but rather 
consider the condition of large values of |£|. Differentiating ( fl4l) one can 
write 

= ^A (17) 
5J(y) 6J(y)^ [ ' } 

We will not write at this point the explicit expression of 8W/SJ(y) but just 
point out its pertinent properties: 

i) It is a linear functional of J and a linear function of £. 

ii) Every derivative of 5WJ SJ(y) with respect to J deletes a factor £. In 
contrast, every derivative of A brings up a factor £. 
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With these considerations we now examine a generic derivative in the sum 
over n in (|T2|) and write in a notation symbolizing the n-order derivative of 
the product in the r.h.s. of ([Tj 



5™ (5A/8J(y)) 



8J(y) 8J(x 1 )...5J(x n ) 



E 

p=0 



^(n-p) 



<5 P 



My) 



(18) 



By property i) the sum over p runs only from to 1, so that we can write 
the expansion (fl2l) as 



5W 



SJ(y) 

n 



J=0 



y2^K( Xl ,...,x n )*A^ 

* — ' n,\ 

K(xt, ...,x n )* 



J=0 



1 



(n-1) 
6J(x n )SJ(y)_ 



By property ii) each "n" term in the first series is two orders in £ higher 
than the corresponding term in the second series and, besides, according to 
expression ffTTT) the first series is equal to P(£). In this way, the sum of 
the second series must be equal to (Ti(y,£) — dW/dJ(y)\ J=0 )P(£), which, 
by construction, is a finite magnitude. We have thus two convergent series, 
one of which is term by term much larger than the other in the limit |£| — > 
oo. Although this is not sufficient to ensure that the sum of the first series 
dominates over the sum of the second series, we obtain a plausible, simple 
model if we assume it and write in the large |£| limit 



r 1 (y,0^8W/6J(y)\ 



A -1 £E(y, xi) * h(xx) 



(19) 



where A = T l {x\,X2) * h{x\)h{x2) ■ The function k does not appear in ([19!) 
because it contributes only to lower order terms in SW/5J(y)\ J=Q . 
The same argument applies to a generic non-linear term as 



r m (yi,-,y m ,€) = {<t>{yi)-<l>{ym)\h{xi) *<f>{x 1 ) =0 

which can be obtained from (compare with Eq. (|12[) ) 



(20) 



r m P(0 



i 

y_]—K{x u ...,x n ) * 



S m A {n) 



5J(y l )...5J{y n 



j=o 
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Properties i) and ii) allow this series to be split into m + 1 series, such that 
m of them are term by term much smaller (by factors f -2 , f~ 2m ) than 
the "dominant" series 



5 m W 



5J(yi)...5J(y r , 



7=0 n 



Again, by comparison with expression (TTTj) . this series is convergent, so that 
the remaining series is also convergent. We thus model T m for large values 
of |f | as 

r ( t\ 5 ™ W 

r m {yi,--,ym,0 - 



5J(yx)...5J(y m ) J=0 

m 5 -Km 

) * h(x m ). (21) 

An important consistency check is that contraction of (|21j) with h{y\)...h(y ri 
satisfies the exact identity 

r m (yi, ...,y m ,0 * h(yx)...h(y m ) = f"\ 

for any h(x\). 

So far no conditions have been imposed on the functions k(x) and a(x, y). 
The only restriction is that both, J r o[4>} and F[(f)\ are pdF's of (f>, so that they 
are normalized to unity for the same measure T><f) 



J F [<f>] Vcf> = J ^[<f>]V<f>=l, 



and that, of course, is independent of k(x) and a(x,y). From these 

conditions and the expression ([6]) of J-"o[(j)] we easily obtain 

W) ~ (22a) 



Sa(x,y) Sa(x,y) 2 

where (...) stands for the average using JF . To choose the functions n(x) 
and a(x,y) consider that we have modeled the large f asymptotics of the 
Tm terms as the sum of a "dominant" series in the sense described above. 
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Since the sums involve infinite series, the term by term dominance is not 
sufficient to ensure dominance of the sum, and that is why what we obtain is 
an approximation whose validity cannot be ascertained in general. However, 
we expect the approximation to be better the faster the series (jTJ) converge, 
so that fewer terms are required in the sums to approximate the functions 
at a given value of £. With only the freedom to choose k(x) and a(x,y) we 
can then expect better results if these functions are chosen so as to keep Q[4>] 
as small as possible in some averaged sense. Taking into account that Q is 
a definite positive magnitude, we could minimize 

or f T\w.Q T>4>, etc.. Of these possible choices, the one that weights pref- 
erentially the small values of Q is j T\nQ T>(j), which, using relations (I22|) . 
leads immediately to ((...) represents the average using T\ that is, the true 
average) 

(cj>{x)) = (<j>(x)) = -E(x, Xl ) * «(x x ), (23) 

(<t>(x)<i>(y)) = (<Kx)<f>(y)) = 

y) + k(zi)k(j/i) * E(ari, x)E(y x , y). (24) 

These equations are easily inverted to determine the functions E = a -1 and 

Kj cLS 

S(x, y) = {mm) - (0(*)> (4>(y)) ■ (25) 
k(x) = -<t(x,xi) * (<f>(xi)) , (26) 
We then summarize the model for the asymptotic averages as 

Fm(yi,-,ym,Q -> C 1 [S(x ,Xq) * h(x )h(x' )f m 

xi) * h(xi)...E(y m , x m ) * h(x m ), (27) 

with T,(x,y) given by Eq. (j2~5]l . 

Finally, the conditional average of the stochastic Gaussian force of zero 
mean can be obtained explicitly for the case of very short lived force corre- 
lation of the forrrP^l 

(f(x,t)f(x',t')) = 5(t-t')F(x-x'), 
and is given in our notation as 

(f(x,t)\0 = - 1 -fh(x 1 )*F(x-x 1 ). (28) 
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III. Applications 



A. Temperature diffusion term 

As a first simple check, we consider the average related to temperature 
diffusion for which measurements exist 8 

s(0 = <V 2 0|0(x) = £) . 

It corresponds to h{x\) = 5 (x\ — x), and can be modeled according to Eq. 
(E2D as 



s(0 = vjr^.o 



y=x 



(0(x)V 2 0)-(0(a;))V 2 
<0 2 (*)) - (<K*)> 2 

It is convenient to take as a non-dimensional field of zero mean and unit 
variance, defined from the physical temperature field T(x) as 

= T(,)-<T(,)) 

where cry is the variance of T. The measured non-dimensional conditional 
average is defined as 

^^((iv^iy 1 ^), (si) 

which, from (1291) and the fact that by definition the variance of is constant, 
reduces to 

r(0 - (32) 

This simple result of the model has been experimentally seen to hold for a 
wide range of values of £P 



B. Forced Burgers turbulence 

For a randomly forced Burgers flow, the velocity u(x,t) satisfies^ 

du du d 2 u „, , j i 

« + "* = "a? + '<*■*>• (33) 
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where v is the kinematic viscosity, f(x,t) is a random force with Gaussian 
distribution of zero mean and variance given by 

</(z, t)f(x\ 0) = - ^ - *'), (34) 

with F an even function of its argument which decays sufficiently fast for 
\x — x'\ larger than a correlation length L c . Besides, -F(O) = 2e', where e' 
is the rate of injection of energy density. Eq. fl33l) with forcing defined 
by (jSH) has been extensively studied using different, powerful theoretical 
method s ' 10 * 11 * 12 * 13 * 14 * 15 ! and numerical simulations] 16 * 17 ^ 



1. Pdf of velocity difference 

The statistical magnitude to be considered first is the instantaneous ve- 
locity difference across a separation r, G[u] = u(x + r,t) — u(x,t). The 
evolution equation for the corresponding pdf P(£, r, t) is, from Eqs. fl2]) and 



dP d d 2 P 
¥ ^((^r)P) + [f(0)-F(r)]^, 



where 



A 



d d 2 
u(v, t)-Q-u(r], t) - v—u{r], t) 



(35) 



(36) 



h{rj) = 8G/8u(rj) = 8(x + r — rf) — 8(x — rj), with (...|^,r), as above, a 
short notation for (...\u(x + r, t) — u(x, t) = £). Finally, statistical spatial 
homogeneity allows to simplify the expression of A to write fl35l) as 



dP 



d_ 



d 



[F(0) - F(r)] 



d 2 u 
dx 2 

d 2 p 



\Z,r)P 



2 ' 



(37) 



Using the model ([27)1 we now evaluate the asymptotic expressions of the 
conditional averages appearing in this equation, which are conveniently ex- 



10 



pressed as 



' d 2 u 



R=0 



R=0 



Using expression (J27jl with m — 1, 2 one obtains asymptotically 



W^)ie,r) 



2S 2 (r 



[&(#)- S 2 (r-i?)] 



(38) 
(39) 

(40) 



(u 2 (R)\^r) 



4Sg(r) 



S 2 2 (J2) 



[S 2 (r - 

2 S 2 (i?)S 2 (r-i?)] 



(41) 



where it was used that in the spatially homogeneous case considered one can 
write Ti(x, x') = (u(x)u(x')) = (u(x)u(x)) — l/2S 2 (x — x'), with S 2 (x — x') = 
({u(x) — u(x')) 2 ) the second order structure function. To evaluate the R 
derivatives appearing in (|38|) and (1391) one needs to know the behavior of S 2 
for small separation. This can be done because at very small scales, less than 
the dissipative scale set up by viscous effects, the velocity field is a smooth 
function of x and so one hafP^«(x + R) = u(x) + Rdu/dx + 0(R 2 ), which 
leads to S 2 {R) = R 2 ((du/dx) 2 ) + 0(R 3 ). With all this, ([38]) and pi} can 
be readily evaluated as 



"2 f 



,d\nS 2 

dr 



/cPu 
\ dx 2 



\C,r 



eU(r) 



(42) 
(43) 



where it was used that v ((d x u) 2 ) = e, with e the rate of dissipation of energy 
density, and where n(r) = 1 — v (d 2 S 2 /dr 2 ) /(2e). The asymptotic equation 
then reads 



dP 
~dt 



1 d\nS 2 d 

2 dr d£ 
2eU d 



(ep) + 

m + [F( o) - F (r)} 



d 2 P 

w 



(44) 
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This equation is not strictly closed because the evaluation of S 2 (r) = J £ 2 P(£, r)d£ 
requires the knowledge of P in the whole range of £, not only in the large £ 
asymptotic where (jUJ) holds. However, 5*2 is a much simpler object to deal 
with, either theoretically or experimentally, than conditional averages. Be- 
sides, much and important information can be gathered without knowledge 
of its explicit expression as shown now. Let us consider now the stationary 
solutions of (|44p . in which case the rate of energy injection e' equals the rate 
of energy dissipation e. The forcing correlation is then written as 

F{r) = 2e X (r/L c ), (45) 

where x is an even, dimensionless function that rapidly decays for arguments 
larger than one, also satisfying x(0) = 1. For r large compared to the forcing 
correlation length L c one then has F(r) ~ 0, and S 2 — 2 (u 2 ), which is 
r independent. With these considerations, (jUJ) immediately leads to the 
Gaussian solution P(£) = exp(— £ 2 / (2S 2 )), as it must since the two velocities 
involved are statistically independent for the r's considered. Much more 
interesting is the situation r « L c , in which case one can write F(r) = 
2e [1 + x"(0) (r/L c ) 2 /2 + 0((r/L c )% with x"(0) < 0. If £ 2 P -> at large 
£, the immediate integral of (jUJ) gives 



T 2 

InP- 



r 2 x"(0) 



ed\ns 2 + £ 2 n 



6e dr S- 



2tt1 

(46) 



If in the range of r considered S 2 (r) behaves as a power law, (in fact, S 2 (r) ~ 
r 2 for r smaller than the dissipation scale, as seen above, and S 2 (r) ~ r in 
the inertial range^) we have d\nS 2 /dr ~ r _1 . Taking into account that 
x"(0) < 0, it then results that the leading behavior for large £ is P ~ 
exp[— L 2 ^ 3 /(er 3 )]. This solution is of course valid only for £ > 0. For £ < 
then, £ 2 P cannot approach zero as £ — » — cxd, and so P must behave as 

P=^, (47, 
where C(r) is an unknown, positive function of r. 

2. Pdf of velocity derivative 

Let us consider now the pdf of velocity derivative at the origin, G [u] = 
du/dx\ x=0 , so that ^(77) = 5G/8u{rj) = —S'(r]), where the prime indicates 
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derivative with respect to the argument. From Eqs. 
obtain 



and (I2"gj) we readily 



dP 
~dt 



d_ 



a 3 



u , 



dx 3 



\t))P 



1 d 2 P 
-i F " (0)^4-. 

2 V ; d£ 2 



We now write, using expression (1271) with m — 2, 

<9 2 



<9ar 



(m(O)m(x) |0 



cc=0 



M — If 



where the zero comes from the assumed statistical homogeneity. Using now 
expression (j27p with m = 1 we have, using also the condition of homogeneity, 



(u—\ l du a2 " X 

2 > dx/ \ dx dx 2 
L ' ' du^ 9 = °' 



<9 3 W 

<9x 3 



<9x 3 



(« (*) 10 



\dx) 



x=0 



so that, calling k = v Ud 2 u/dx 2 f){(du/dxf) 



yd 2 u/dx 



2\ 2 



we 



have for the statistically stationary case, using also the expression (H5|) for 
the forcing, 



d_ 



x "(0) OP 



0. 



LI di 

and analogously to the velocity difference case we obtain for £ > 



and, for f < 0, 



X"(0)V3£ 2e 



All these results for the Burgers equation reproduce well the behavior 
obtained theoretically and numerically in the references cited above. 
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IV. Conclusions 



We have presented a simple closure for the asymptotic pdf evolution equa- 
tion of rather generic turbulence problems. We have presented results for the 
simplest possible problems in turbulence, for which a rather detailed knowl- 
edge exists. For the paradigmatic problem of Navier-Stokes turbulence the 
application of the formalism is more involved, and has been applied in a par- 
ticular version inP where very reasonable results were obtained, although 
comparison with detailed pdf 's is more difficult. 
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